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We study the statistics of pair dispersion in two-dimensional turbulence. Direct 
numerical simulations show that the pdf of pair separations is in agreement with 
the Richardson prediction. The pdf of doubling times follows dimensional scaling 
and shows a long tail which is the signature of close approaches between particles 
initially seeded with a large separation. This phenomenon is related to the formation 
of fronts in passive scalar advection. 
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1 Introduction 

The concentration of a dilute solution of a passive tracer in an incompressible 
flow obeys the scalar equation 



where v is the divergenceless velocity field, k is the molecular diffusivity, and 
/ is the external source of tracer fluctuations. Equation (1) can be solved by 
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the method of characteristics to obtain the solution 

0(x,f) = f^ds f(p(s),s) ^ 
p(s) =v(p(s),s) + V2kti(s), p(t)=x. 

The characteristics p(t) are the Lagrangian trajectories of the fluid particles, 
and y/2nr] is the white noise contribution due to molecular diffusion. There 
is an immense literature devoted to both Eulerian, (1), and Lagrangian, (2), 
descriptions of passive transport in turbulence [1]. What is relevant to our 
purposes is to keep in mind the tight relationship between these two comple- 
mentary descriptions. As an instance, simultaneous two-point correlations of 
the scalar field can be written in terms of two-particle Lagrangian statistics 
as 

t t 

(0( Xl ,i)0(x 2 ,t))= J d Sl J ds 2 (f(p(s l ),s 1 )f(p(s 2 ),s 2 )). (3) 

— oo — oo 

Properly choosing the form of the correlation function of the scalar forcing, e.g. 
(/(xi, £i)/(x 2 , t 2 )) = x(|x x — x 2 |)<5'(ti — ^2)5 and exploiting space homogeneity, 
expression (3) can be further simplified to the form 

t 

(0(x,t)0(x + R,t)) = J ds J dr X (r)p(r,s\K,t) . (4) 

—00 

where p(r, s|R, t) is the probability density function for a pair to be at a 
separation r at time s, under the condition that it has to have a separation R 
at time t. 

Summarizing, the knowledge of the statistics of pair dispersion is sufficient to 
determine the values of the correlations of passive scalar for a given forcing. 
On the contrary, to extract the Lagrangian statistics from the Eulerian one it 
is necessary to know the scalar correlation functions for different forcings. In 
this sense, the Lagrangian information is more fundamental, and hereafter we 
shall restrict to this one. 



2 Statistics of pair dispersion 

The dispersion of a particles' pair in turbulence can be phenomenologically 
described in terms of a diffusion equation for the probability distribution of 
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pair separations 



dpjr, t) 
dt 




K(r,t) 



dyij, t) 
drj 



) 



(5) 



with a space and time dependent diffusion coefficient K(r,t). In general, the 
description by means of a diffusion equation is a drastic simplification. The 
only case in which it can be proven that particle separations obey to (5) is 
when the advecting velocity field is rapidly changing in time [2,3]. 

The original Richardson proposal, obtained from experimental data in the 
atmosphere, is K(r,t) = K{r) ~ r 4 / 3 [4]. This leads to the well known non- 
Gaussian distribution 



From the Richardson distribution (6) one has immediately that the mean 
square particle separation grows as 



The "t 3 " law, which is known as the Richardson law, has been observed, 
although with large uncertainty, in direct numerical simulations [5] and, more 
recently, in laboratory experiments [6]. 

The diffusion equation (6) is not the unique possibility which leads to the 
"t 3 " law. Batchelor [7] assuming that the diffusion coefficient should depend 
on average quantities, proposed that K(r,t) = K(t) ~ (r 2 (t)) 2//3 ~ t 2 . In 
this case the distribution of pair distances is Gaussian with a superballistic 
growing variance 



Of course, this is not the end of the story. Formally, any diffusion coefficient 
of the form K(r, t) ~ r a t b with 3a + 2b = 4 is compatible with the "t 3 " law 
but gives different distribution function [1,8]. Early experimental data were 
in favor of the Batchelor Gaussian distribution (8) [1], but recent laboratory 
experiments are more in the direction of the Richardson original proposal [6]. 

The Richardson law can be derived by a simple dimensional argument which 
makes use of the Kolmogorov similarity law for the Eulerian velocity incre- 
ments in fully developed turbulence. By definition one has 




(6) 



R 2 (t) = (r 2 (t)) ~ t z 



(7) 




(8) 



d 1 
~dt2 



R 2 (t) = (r-oV L )(r)> = (r5v^\r)} 



(9) 
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where <5v( L )(r) = v(x(£) + r) — v(x(£)) is the Lagrangian velocity increment 
and 5v\\ is its projection on the r direction. Assuming Kolmogorov scaling 
for Lagrangian velocity differences, 6v\\(r) ~ r 1 / 3 , one obtains from (9) the 
Richardson law (7). The assumption that the Lagrangian velocity difference 
has the same Kolmogorov scaling of the Eulerian one relies on the intuitive 
idea that the main contribution to the separation rate follows from eddies with 
size comparable with the separation itself. 

Let us conclude this section with a remark. The description of relative disper- 
sion in terms of the diffusion equation (5) assumes a self-similarity in the pro- 
cess. Of course this could not be the case, e.g. in intermittent three dimensional 
fully developed turbulence. As a matter of fact, in presence of intermittency 
of the velocity field, one can expect a kind of "Lagrangian intermittency" , in 
the sense that different moments (r p (t)) have different scaling exponents: 

(r p (t)) ~ t°* (fO) 

with a p 3/2p. This problem has been discussed in several papers [9-12] with 
different conclusions. Recent detailed investigations with a synthetic turbulent 
model gave the evidence of Lagrangian intermittency with scaling exponents 
a p linked to the Eulerian intermittent scaling exponents [13]. 



3 Pair dispersion in two-dimensional turbulence 

Pair dispersion statistics has been investigated by direct numerical simulation 
of the inverse energy cascade in two-dimensional turbulence. There are several 
reasons for considering 2D turbulence. From an applicative point of view, 
two-dimensional Navier-Stokes equations are among the simplest systems of 
geophysical interest. The observed absence of intermittency [14,15] makes the 
2D inverse energy cascade an ideal framework for the study of Richardson 
scaling. Moreover, the dimensionality of the problem makes feasible direct 
numerical simulations at high Reynolds numbers. 

The 2D Navier-Stokes equations for the vorticity uo = V x v = — Aip are: 

d t ui + J (uj, ip) = uAu — auj + 0, (11) 

where ip is the stream function and J denotes the Jacobian. The friction linear 
term -aw extracts energy from the system to avoid Bose-Einstein condensa- 
tion at the gravest modes [16]. The forcing is active only on a typical scale If 
and is ^-correlated in time to ensure the control of the energy injection rate. 
The viscous term has the role of removing enstrophy at scales smaller than 
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Fig. 1. Energy spectrum E(k). The dashed line is the Kolmogorov scaling 
E(k) ~ /c -5 / 3 . In the inset it is shown the energy flux H(k). 

If and, as customary, it is numerically more convenient to substitute it by a 
hyperviscous term (of order eight in our simulations). Numerical integration 
of (11) is performed by a standard pseudospectral method on a doubly peri- 
odic square domain of N 2 = 2048 2 grid points. All the results presented are 
obtained in conditions of stationary turbulence. 

In Figure 1 we present the energy spectrum, which displays Kolmogorov scal- 
ing E(k) = Ce 2//3 /c" 5//3 over about two decades with Kolmogorov constant 
C ~ 6.0. The inertial range correspond to the region of constant flux, also 
plotted in Figure 1. Previous numerical investigation has shown that veloc- 
ity differences statistics in the inverse cascade is almost Gaussian with Kol- 
mogorov scaling not affected by intermittency corrections [15]. In this case 
we expect the Lagrangian statistics to be self-similar with Richardson scaling 
[13]. 

In Figure 2 we plot relative dispersion R 2 (t) obtained after averaging over 
64000 particle pairs for two different initial conditions R 2 (0). The Richardson 
t 3 law is observed in a limited time interval, especially for the largest R 2 (0) 
run. It is remarkable that the relative separation law displays such a strong 
dependence on the initial conditions even in this high resolution runs. 

The probability density functions of pair separations is plotted in Figure 3 at 
two different times. For short times (t = 0.2) in which the relative dispersion 
is in the Richardson regime (see Figure 2) we see that the Richardson distri- 
bution (6) fits well the numerical data. This is, we think, a clear evidence of 
the substantial validity of the original Richardson description. Let us observe 
that until now this point was not clear: recent laboratory data [6] pointed 
to the Richardson distribution but there were strong deviations from (6). At 
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Fig. 2. Relative dispersion R 2 (t) for two initial separation R(0) = 1.5 x 10 3 (+) 
and R(0) = 3 x 1CT 3 (x). The continuous line is the Richardson law R 2 (t) ~ t 3 . 
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Fig. 3. Probability density function of relative separation at times t = 0.2 and 
t = 5.0 rescaled with R(t) = (r 2 (t)) 1 / 2 . The continuous line is the Richardson 
prediction (6), the dashed line is a Gaussian distribution. 

long time t = 5.0 relative separation distribution is described by a Gaussian 
distribution, but this has no relation with the Batchelor proposal (8) because 
it is not in the scaling range. At time t — 5 the average separation is of the 
size of the computational box and we have normal dispersion a la Taylor. 

The Richardson pdf has a strong cusp at r = which signals the high prob- 
ability for a pair to reach a very small separation compared to the typical 
value R(t). As we shall see later on, this effect is highlighted by considering 
the statistics of first exit-times. 
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4 Doubling time statistics 



We have seen in the previous section that the Richardson picture seems to 
be confirmed in the inverse energy cascade in two dimensions. Nevertheless 
we have seen that it is difficult to observe the Richardson scaling law even 
in our high resolution direct numerical simulations. To understand this effect 
consider a series of particle pair dispersion experiments, in which a couple of 
particles is released at time t — with initial separation R(0). At a fixed time 
t one performs an average over all different experiments and computes R 2 (t). 
It is clear that, unless t is large enough that all particle pairs have "forgotten" 
their initial conditions, the average will be biased. This is at the origin of the 
flattening of R 2 (t) for small times, which we can call a crossover from initial 
condition to self similar regime. Of course, at larger R(0) correspond longer 
crossover regimes (see Figure 2). A similar effect is observed for times of the 
order of the integral time-scale since some particle pairs might have reached 
a separation larger than the integral scale and thus diffuse normally, biasing 
the average, so that the curve R 2 (t) flattens again. 

To overcome this difficulty we use an alternative approach based on statistics 
at fixed scale [17]. The method has been successfully applied to the analysis of 
relative dispersion in synthetic turbulent flow [13] and in experimental convec- 
tive laminar flow [18]. The method works as follows. Given a set of thresholds 
R n = r n R(0) within the inertial range, one computes the "doubling time" 
T r (R n ) defined as the time it takes for the particle pair separation to grow 
from threshold R n to the next one R n +i- Averages are then performed over 
many dispersion experiments, i.e., particle pairs. The outstanding advantage 
of this kind of averaging at fixed scale separation, as opposite to a fixed time, 
is that it removes crossover effects since all sampled particle pairs belong to 
the inertial range. 

The scaling property of the doubling time statistics in fully developed turbu- 
lence is obtained by a simple dimensional argument. The time it takes for the 
particle pair separation to grow from R to rR can be dimensionally estimated 
as T r (R) ~ R/5v^ L \R); we thus expect for the inverse doubling times the 
scaling (T r (R)} ~ i? 2//3 . From the definition, doubling times depend on the 
threshold ratio r. It is then useful to consider the normalized quantity 



m= wm loer (12) 



which is called the Finite Size Lyapunov Exponent (FSLE) because is reduces 
to the (Lagrangian) Lyapunov exponent in the limit of small separation R — > 
[19]. 
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Fig. 4. Finite Size Lyapunov Exponent (12) for the same trajectories of Figure 
2. The initial threshold is R(0) = 0.0031 and the ratio is r = 1.2. The line is the 
theoretical Richardson scaling i?~ 2 / 3 . 

In Figure 4 it is shown the FSLE for the same simulation of Figure 2. At 
small scales R < 0.01 there appears a constant plateau corresponding to the 
Lagrangian Lyapunov exponent X(R) — 19. At larger R, we observe the power 
law X(R) — R~ 2 ^ 3 on a scaling range which is well enhanced with respect to 
the relative dispersion of Figure 2. 

The scaling of doubling times gives also information about the two-point cor- 
relations of passive scalar. Indeed, assuming that the correlation of the scalar 
forcing decays rapidly to zero with a typical scale L, by means of equation (4) 
we obtain 

t 

(0(x,f)0(x + R,f)) ~ J ds j drp(r,s|R,t) (13) 

-oo |r|<L 

that is the average time that a particle pair released at a separation R spends 
below the scale L. Thus, our results for doubling times enable us to express 
the scalar correlation as 

(0(x, t)9(x + R, t)) ~ L 2 ^ — R 2 / 3 (14) 

for any L and R in the scaling range of the velocity field. Direct numerical 
simulations of passive scalar advection, eq. (1), in the Navier-Stokes flow gen- 
erated by eq. (11) confirm that the exponent of scalar correlations is indeed 
indistinguishable from 2/3 [20]. 



Beyond scaling properties of averaged quantities, the inspection of the pdf of 
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Fig. 5. Probability distribution functions of doubling times rescaled with the average 
doubling time (T(R)) for different R in the inertial range. 

doubling times is very insightful to capture the main features of pair dispersion. 

In figure 5 it is shown the pdf of doubling times rescaled with their aver- 
age values. The normalized pdf's at different separations in the inertial range 
collapse, indicating the self-similarity of the Lagrangian dispersion. Most im- 
portant, there is a large number of events for which the pair wanders for 
20 — 30 times the average value before exiting. This effect is a reflection of 
the strong cusp observed in the pdf of pair separations. In these events the 
particles which are initially at a separation lying well inside the inertial range 
can approach each other as close as the diffusive scale. In the language of 
the passive scalar field, since trajectories originating from widely separated 
regions of space can carry very different values of the concentration field 9, 
these approaches generate steep gradients of scalar across small scales. These 
structures, known as "cliffs" , have been actually observed both experimentally 
[21-23] and numerically [24,25] for the temperature field in a turbulent flow. 



5 Conclusion 



Passive scalar transport in turbulence can be described in two complementary 
ways. It is possible to adopt the field description, and think in terms of cor- 
relation functions of the scalar field, either to prefer the particle description, 
and thus ask questions about the statistics of relative separations. These two 
aspects complete each other. We have investigated the Lagrangian properties 
of transport in two-dimensional turbulence, our results show that pair disper- 
sion statistics is not intermittent, since the velocity field is self-similar and the 
geometric content of two-particle configurations is trivial. This result makes 
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contact with the analytical result that can be derived in the rapid-change 
model, where two-point scalar correlations show no anomaly. This will not be 
the case in three-dimensional turbulence, a case which is under current anal- 
ysis. 

Although on average particle pairs separate, probability density functions of 
pair separations and of doubling times clearly display the fingerprint of fre- 
quent close approaches between particles. These events occur with a relatively 
high probability, and are for the formation of quasi- discontinuities in the scalar 
field (the cliffs). 
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